Multiple and virtual photon processes in radiation-induced magnetoresistance 
oscillations in two-dimensional electron systems 
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Recently discovered new structures and zero-resistance states outside the well-known oscillations 
are demonstrated to arise from multiphoton assisted processes, by a detailed analysis of microwave 
photoresistance in two-dimensional electron systems under enhanced radiation. The concomitant 
resistance dropping and the peak narrowing observed in the experiments are also reproduced. We 
show that the radiation-induced suppression of average resistance comes from virtual photon effect 
and exists throughout the whole magnetic field range. 



The recent interest in radiation related magneto- 
transport in two-dimensional (2D) electron gas 
(EG) has been stimulated by the discovery of mi- 
crowave induced magnetorersistance oscillations 
(MIMOs) and zero-resistance states (ZRS) in ultra- 
high mobility systems at low temperatures^* 2 * 3 ^ 
Tremendous experimenta l^?i'i'i?i?i 1 9i 11 i 12 i 1 ?] 14 and 
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efforts have been devoted to study this exciting phe- 
nomenon and a general understanding of it has been 
reached. MIMOs emerge as the magnetoresistance R xx 
of the 2D system subject to a microwave radiation of 
frequency / = uj/2it, exhibiting periodic oscillation as a 
function of the inverse magnetic field l/B. It features 
the periodical appearance of peak-valley pairs around 
lu/lo c — 1, 2, 3, 4, • • •, i.e. a maximum at lo/lu c = j — 5j 

and a minimum at uj/oj c = j + Sj , with j = 1, 2, 3,4, • • • 
and < 5^ < 1/4. Here u> c is the cyclotron frequency 
and lu/uj c = j are the node points of the oscillation. 
With increasing the radiation intensity the minimum 
value of R X x drops down until a vanishing resistance is 
measured, i.e. ZRS. 

In addition to these well-established features, a clear 
secondary structure between the first and second main 
peaks was observed experimentally* 3 "*^ and predicted 
theoretically 3 , with reference to two-photon processings 3 . 
A distinct minimum and/or maximum outside the first 
main peak was also observed experimentally at 47.5, 52 
and 58GRz^iL at 45 and 35GH3*U, and at 30GH3& 
which was referred to two-photoniAi or the second har- 
monic effect^ Theoretically a similar structure due to 
two-photon process was anticipated at 60 and 40 GHz, 23 
followed by further experimental observations at lower 
frequency (10-50 GHz). SiiLH A more prominent exper- 
imental result was recently reported by Zudov et al at 
27 GHzfiA where, with intensified microwave radiation 
the above-mentioned minima evolve into ZRS. Despite 
the consensus on the existence of these structures the 
assignment of their locations has so far been different 
among different groupsiii Physically these peak-valley 
structures can be due to multi-photon procesSfiiiiMiSiSi 
or due to higher harmonics^ of the radiation-induced 
high-frequency current. To ascertain these structures as 
arising from multiphoton effect and to help identifying 



their positions, it is imperative to reproduce them con- 
vincingly from a careful theoretical calculation extended 
to higher radiation intensity. 

Another feature is the descent of the average dis- 
sipative resistance under microwave irradiation, which 
was experimentally observed on the high magnetic field 
side where MIMO shows up relatively weak&i^i^i^ 
Theoretically, the magnetoresistance drop was also 
anticipated*^* 3 ^ and referred to virtual photon effect^ 
or to microwave-induced dynamic localization* 3 ^ Never- 
theless, the range and the degree of this resistance de- 
scent remain unclear, and a unified theoretical treatment 
covering higher radiation intensity to demonstrate it to- 
gether with multiphoton structures is desirable. 

In this Letter we report on a detailed analysis of MI- 
MOs with enhanced radiation intensity, based on a the- 
oretical model which covers all orders of real and vir- 
tual photon processes of the base frequency, but excludes 
higher harmonic current response. 

The model considers that a dc electric field Eo and a 
high frequency (HF) field E(£) = E s sin(wi) + E c cos(ujt) 
are applied in a quasi-2D system consisting of N e in- 
teracting electrons in a unit area of the x-y plane, to- 
gether with a magnetic field B = (0, 0, B) along the z 
direction. The approach is based on the separation of 
the center-of-mass motion from the relative electron mo- 
tion of the electrons and describes the transport state 
of a high-carrier-density many-electron system under a 
radiation field in terms of a time-dependent electron 
drift velocity oscillating with base radiation frequency, 
v(i) = v c cos(wt) + v s sin(wi), and another part v , de- 
scribing the slowly varying electron drift motion, together 
with an electron temperature T c characterizing the elec- 
tron heatingi* 3 ^*^ In the case of ultra-clean electron gas at 
low temperatures, Vo and T e satisfy the following force- 
and energy-balance equationsi^iill 

m^= e E + e (v xB) + ^, (1) 

N C E • v + S p - W = 0, (2) 

with v c and v s determined by 

-mwv c = eE s + e(v s x B), (3) 
wv s = eE c + e(v c xB). (4) 
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Here e and m are the electron charge and effective mass, 

oo 

F =]T|[/(q||)| 2 £ q|,^(0n 3 (q|,,wo-na;) (5) 

q ii n= — oo 

is the damping force of the moving center of mass, 

oc 

5 p = ^|C/(q||)| 2 nwJ2(0II 2 (q|| )W d-nw) ( 6 ) 

q|| n— — oo 

is the averaged rate of the electron energy absorption 
from the HF field, and 

oo 

TY = ]T|M(q)| 2 « q J 2 (OA 2 (q^ + r! q -^) (7) 

q n— — oo 

is the average rate of the electron energy loss to the 
lattice. In the above equations, J n (0 is the Bessel 
function of order n, £ = v/(qjj ' Vc ) 2 + fall ' Vs ) 2 / W ! 
ci)q = q|| • vq, t/(q||) and M(q) are effective impurity and 
phonon scattering potentials, lis (tjll 5 f2) and A 2 (q, f2) = 
2n 2 (q||,0)[n(fi q /T)-n(rJ/T e )] (with n(x) = l/(e*-l)) 
are the imaginary parts of the electron density correla- 
tion function and electron-phonon correlation function of 
the system in the magnetic field. 

The IT2 (ct 1 1 function can be expressed in the Landau 
representation: 37 

n 2 (q||,^) - ^£C n , n ,(Z|gjf/2)II 2 (rc ) n , ,fi), (8) 

^ n,n' 

n 2 (n, n', n) = - ~ /" <fe [/(e) - /(e + fi)] 

xImG n ( £ + 0)ImG n .(£), (9) 

where C n , n+J (Y) = n![(n + i)!]- 1 ^-*^^) ] 2 with 
(Y) the associate Laguerre polynomial, Ib = y/TJ\eB\, 
/(e) = {exp[(e — /tx)/T e ] + is the Fermi function at 
electron temperature T c . The density of states of the 
n-th Landau level is modeled with a Gaussian form»2& 

ImG„(e) = -(v^F/T) exp[-2( £ - e„) 2 /F 2 ]. (10) 

having a half- width T = (SeLUca/Trmno) 1 ^ 2 around the 
level center e n . Here no is the linear mobility at lattice 
temperature T in the absence of magnetic field, uj c — 
eB/m and a is a semiempirical broadening parameter. 

For time-independent vo , we obtain the transverse and 
longitudinal dc resistivities from Eq. R xy = B/N e e, 

R xx = -F -v /(N 2 e 2 v 2 ). (11) 

The (linear) magnetoresistivity is its vo — > limit. 
Within certain field range, R xx can be negative at 
small vq, but increases with increasing vq and passes 
through zero at a finite uor& implying that the time- 
independent small-current solution is unstable and a spa- 
tially nonuniform 1 ^ or a time-dependent solution 34 may 




FIG. 1: The magnetoresistivity R xx of a GaAs-based 2DEG 
with N c = 3.0 x 10 15 m~ 2 , fj, = 2000m 2 /Vs and a = 10, 
subjected to 50 GHz radiations of incident amplitudes E- ls = 
4, 5, 6.5 and 8.5 V/cm at lattice temperature T = IK. 



develop, which exhibits measured zero resistance. There- 
fore we identify the region where a negative dissipative 
magnetoresistance develops as that of the ZRS. 

Assume that the 2DEG is contained in a thin sam- 
ple suspended in vacuum at plane z = 0. When an 
electromagnetic wave illuminates the plane perpendicu- 
larly with the incident electric field E; (t) = Ei s sin(wt) + 
E ic cos(u;i), the HF electric field in the 2DEG is 

E* = % U +Eift). (12) 

Using this E(t) in Eqs. © and Q, vi and v 2 are explic- 
itly expressed in terms of incident field E; s and E; c . 

Note that all orders of real (\n\ > 0) and virtual (n = 0) 
photon processes are included in the summations over n 
in Eqs. JSJ, © and J7J). The phonon contributions to F 
and S p have been neglected because of the low tempera- 
ture setup in the experiments. The short-range scatterers 
are considered to give the dominant contribution to the 
resistance and energy absorption for the ultra-clean sam- 
ples used^i The numerical calculations are performed for 
x-dircction (parallel to Eo) linearly polarized incident mi- 
crowave fields [Ej S = (Ei, 0), Ei C = 0], using the material 
parameters of GaAs^ 

Figure 1 shows the calculated magnetoresistivity R xx 
as a function of lo/lo c for a GaAs-based 2D system hav- 
ing electron density iV c = 3.0 x 10 15 m~ 2 , linear mobility 
/io = 2000m 2 /Vs and broadening parameter a = 10, 
irradiated by 50 GHz microwaves of three incident am- 
plitudes Ei = 4,5 and 6.5 V/cm at lattice temperature 
T = 1 K. We see that in addition to the main maximum- 
minimum pairs at uj/lu c — 1,2,3 and 4, prominent peak- 
valley pairs show up around uj/lu c — 1/2,2/3,3/5,2/5 
and 2/7, and the minima of pairs lu/lo c = 1/2 and 
uj/u> c = 3/5 drop down to negative. 

Careful analysis shows that the appearance of oscilla- 
tory magnetoresistance comes from real photon-assisted 
electron transitions between different Landau levels as 
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FIG. 2: The Magnetoresistivity R xx of a GaAs-based 2DEG 
with No = 3.0 x 10 15 m~ 2 , jito = 2000m 2 /Vs and a = 5, 
subjected to 27 GHz radiations of incident amplitudes Ei s — 
1.5, 2, 3 and 4 V/cm at lattice temperature T = 1 K. 



indicated in the summation of the electron density- 
correlation function in Eq. JSJ. We denote a real-photon 
assisted process in which an electron jumps across I 
Landau-level spacings with the assistance (emission or 
absorption) of n photons as nw:Zw c , or n:l. This process 
contributes, in the R xx vs u> c /u> curve, a pair structure 
consisting of a minimum and a maximum on both sides 
of lUc/lu = n/l. The location of its minimum or maxi- 
mum may change somewhat depending on the strength 
of the incident microwave, but the node point, which 
is roughly in the center, essentially keeps at the posi- 
tion lo c /lo = n/l. Therefore we use its node position, 
rather than its minimum or the maximum, to identify 



a pair structure. Thus, the single-phonon process 1:1, 
the two-phonon process 2:2, the three-photon process 
3:3,- • •, all contribute to the minimum-maximum pair 
around uj c /uj = 1; the single-phonon process 1:2, the 
two-phonon process 2:4, the three-photon process 3:6,- • •, 
all contribute to the minimum-maximum pair around 
u) c /us = 1/2; etc. These are indicated in Fig. 1, as well as 
in Fig. 2, where we plot the evolution of magnetoresistiv- 
ity R xx with increasing radiation intensity E\ = 1.5,2,3 
and 4 V/cm vs u> c /ui, for a GaAs-based system with 
N c = 3.0 x 10 15 m~ 2 , fi = 2000m 2 /Vs and a = 5, irra- 
diated by 27 GHz microwaves. Several other valley-peak 
pairs associated with 2-, 3-, and 4-photon assisted pro- 
cesses are also identified. 

The predicted two-photon structure 2:1 centered at 
lu c /uj = 2 with a minimum around uj c /lu ~ 1.73 (Fig. 2(a) 
and (b)) reasonably compares with the second harmonic 
minimum around u> c /u> = 8/5 observed in Refs.00- The 
two-photon structures 2:3 centered at uj c /uj = 2/3 and 
2:5 centered at <jj c /u = 2/5 are in agreement with exper- 
imental findings A^iii 

Note that, with the progressive emergence of new 
multiphoton-related pairs when increasing radiation 
power, the peaks (valleys) of the low-order photon re- 
lated pairs become narrower, as is clearly seen for the 
single-photon related pairs at uj/lu c = 1,2 and 3 in both 
figures. This feature and the anticipated positions of the 
peaks and possible ZRSs, agree with the recent experi- 
mental observation^ 

Another interesting aspect is that, concomitantly with 
enhanced R xx oscillation, the average magnetoresistance 
descends down significantly with increasing radiation 
power. This resistance drop is due to effect of virtual 
photon processes, i.e. intra-Landau-level electron scat- 
tering by impurities with simultaneous emission and ab- 
sorption of an arbitrary number of photons. To show this 
we plot the resistivity contributed from the virtual pho- 
ton processes alone, i.e. the n = term (J ) in Eq. JHJ), in 
both figures. The resistance suppression appears almost 
throughout the whole magnetic field range. It may not 
be previously noticed in the region where exhibits strong 
R xx oscillation and ZRS. 
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